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• • . 1 INTRODUCTION 



The weak lensing shear signal has been measured numerically in iV-body simula- 
tions at 14 different redshifts (z s = 0.1 to 3.6) and on angular scales of 9 = 2' to 32'. 
In addition, the data have been validated by analytical computations for an identical 
cosmology, with density parameter Q, m = 0.3 and vacuum energy density parameter 
A = 0.7. This paper reports on the scale and redshift dependence of the shear variance, 

(7 2 ), which may be described by a simple formula of the form {j 2 }(9, z s ) — a(9)z£ S \ 
The redshift dependence for source redshifts up to 1.6, is found to be close to z 2 , which 
is a stronger dependence than earlier analytical predictions (X7 2 ) oc zj- 52 ), although, 
at higher redshifts, the z s dependence of the shear variance is clearly less steep. The 
strong redshift dependence further emphasises the need to know the precise redshift 
distribution for the galaxy sources in any given survey, so that they can be allocated 
to redshift bins accordingly, and the cosmic shear signal correctly interpreted. Equa- 
tions are also given for the variance in the reduced shear, which is a more directly 
measurable quantity observationally. 
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Methods: 



Studies of weak gravitational lensing in cosmology are a very 
powerful tool for attempts to understand the distribution of 
mass and the evolution of the large-scale structure in the 
universe. Recently, it has also become possible to constrain 
values for the cosmological parameters from weak lensing 
studies. Since the gravitational deflections of light arise from 
variations in the gravitational potential along the light path, 
the deflections result from the underlying distribution of 
mass which is usually considered to be in the form of dark 
matter. The lensing signal therefore contains information 
about the clustering of mass along the line-of-sight which 
may be different from the clustering inferred from galaxy 
surveys which trace the luminous matter only. In addition, 
by studying the way light from large numbers of sources at 
high redshifts is deflected, it is possible to obtain informa- 
tion about the way the clustering of mass evolves with time. 

As a result of weak gravitational lensing, a source at 
high redshift will appear magnified (or de-magnified) as the 
beam converges (or diverges) due to matter (or an under- 
density) contained within it. The image also undergoes 
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shearing due to deflections from matter outside the beam, 
and causes a circular source, for example, to appear as an 
elliptical image. Sources at similar redshifts and contained 
within a small field-of-view will display similar magnifica- 
tion and shear characteristics because their light will have 
passed along similar density paths. For this reason there 
will be strong correlations in the changes to the ellipticities, 
particularly on small scales, and declining correlations on 
increasing angular scales. 

The magnitude of this correlation depends strongly on 
the density parameter and the value of the cosmological con- 
stant for the universe, as these parameters reflect both the 
amount of mass and the rate of evolution of structure. A 
number of attempts have been made to estimate weak lens- 
ing statistics in different cosmological models, both analyt- 
ically and numerically. See Bernardeau, Van Waerbeke & 
Mellier (1997), Jain & Seljak (1997) and Bacon, Refregier 
& Ellis (2000), for example, for analytical considerations, 
and Barber, Thomas & Couchman (1999), Jain, Seljak & 
White (2000), Barber et al. (2000), Hamana, Colombi & 
Mellier (2000), Van Waerbeke et al. (2001a) and Premadi et 
al. (2001), for example, for work done using various cosmo- 
logical A-body simulations. These authors have attempted 
to predict values for various lensing statistics, many of which 
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may be obtained directly or indirectly from observational 
data. 

A number of observational results have now been re- 
ported for the so-called cosmic shear signal; see, for example, 
Bacon et al. (2000), Kaiser, Wilson & Luppino (2000), Maoli 
et al. (2001), Van Waerbeke et al. (2000a, b), Wittman et 
al. (2000), Mellier et al. (2001), Rhodes, Refregier & Groth 
(2001) and Van Waerbeke et al. (2001b). There is also a 
number of ongoing and planned observational programmes 
such as the Sloan Digital Sky Survey and surveys using 
the Visible and Infra-red Survey Telescope for Astronomy 
(VISTA). These and other telescope programmes will help 
to further constrain the cosmic shear values and the cos- 
mological parameters. An excellent review of weak gravi- 
tational lcnsing and the measurement of the cosmic shear 
signal in particular is contained in the paper by Bartelmann 
& Schneider (2001). 

The purpose of this paper is to predict important cos- 
mic shear statistics using cosmological TV-body simulations, 
which can be compared directly with observations. In par- 
ticular, it reports the redshift and angular scale dependence 
for the shear variance and the reduced shear, computed nu- 
merically and supported by analytical computations. The 
detailed results are presented for background sources at 14 
different redshifts and angular scales from 2' to 32'. The find- 
ing that the shear variance has a strong redshift dependence 
draws attention to the need to have good redshift data in 
observational surveys before the cosmic shear signal can be 
correctly interpreted in terms of the underlying cosmology. 

The numerical method uses the algorithm described 
in Couchman, Barber & Thomas (1999), which computes 
the three-dimensional shear in the simulations, and which 
has been applied to a simulation with density parameter 
Q. m — 0.3 and vacuum energy density parameter Ao = 0.7. 
(Cosmologies of this type will be referred to as LCDM cos- 
mologies.) To obtain the required statistics on different an- 
gular scales, the computed shear values have been combined 
(using the appropriate angular diameter distance factors 
and accounting for multiple deflections) along lines of sight 
arranged radially from the observer's position at redshift 
z — 0. To support the numerical results, analytical com- 
putations have also been made for an identical cosmology, 
using the procedure of Van Waerbeke et al. (2001a). 

A brief outline of this paper is as follows. 

In Section 2, the equations and definitions required for 
weak lensing are presented, and analytical work to determine 
weak lensing statistics are described. 

In Section 3, the shear algorithm and the V-body sim- 
ulations are summarised together with details of the appli- 
cation of the code for the generation of the lensing statistics 
at the observer's location. 

Section 4 presents the numerical results for the vari- 
ances in the shear and the reduced shear for sources at 
different redshifts, and on different angular scales. Simple 
equations are presented which fit the data to describe these 
variances as functions of redshift and angular scale. 

Section 5 is a discussion of the results. 



2 WEAK LENSING THEORY 

An explicit expression for each element of the three- 
dimensional shear at an arbitrary position, R, arising from 
matter at positions R' has been given by Barber et al. 
(2000): 
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In this expression, (j> is the peculiar gravitational potential, 
Xi and Xj are position coordinates, with the suffixes i and 
j representing the directions denoted by 1, 2 or 3, G is the 
universal gravitational constant, p and p are the density and 
mean density respectively, and a is the scale factor for the 
universe. The three-dimensional shear values computed by 
the shear algorithm can be identified with this expression. 

To compute the required physical properties at z = 
arising from the light of a distant source sheared by a sin- 
gle deflector, the quantities needed are the two-dimensional 
second derivatives of the effective lensing potential, which 
are denoted by ipij , and which are calculated by integrating 
the three-dimensional shear along the coordinate direction 
x;i , and including the appropriate angular diameter distance 
factors: 



DdDds 2_ r d 2 <j>(x 3 ) 
D a c 2 / dxidx 



-dxz- 



(2) 



Here Dd, Dds, and D s are the angular diameter distances 
from the observer to the lens, the lens to the source, and the 
observer to the source, respectively, and c is the velocity of 
light in vacuo. 

Where there are multiple deflections, these quantities, 
evaluated for the ith deflector (i is now used as the deflector 
index) , form the elements of the shear tensor, , equiva- 
lent to the derivatives of the reduced deflection angle, 



U 



■ _ / V>U V>12 
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(3) 



By combining the for all the deflectors, the final prop- 
erties at z = can be evaluated. To do this use is made 
of the multiple lens-plane theory, which has been described 
concisely by Schneider, Ehlers & Falco (1992). In this the- 
ory, Jacobian matrices, which describe how small changes 
in the position vector of an element of the source relate to 
small changes in the position vector of the corresponding el- 
ement in the image, are constructed by recursion. The final 
Jacobian matrix, A, at z = resulting from N deflectors is 
given by 



A = T- 



N 



(4) 



where T is the identity matrix, and in which the individual 
Jacobian matrices are 



A u) = 1 - /3yW (i U (<) 



(5) 



for the jth lens, and 
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for the first lens. In equation |H[ 

= D ± Di i 
Pl] ~ I),- Dj ' 



(6) 



(7) 



where D s and Dj are the angular diameter distances to the 
source and the jth lens respectively, and Di S and Dij are the 
angular diameter distances from the ith lens to the source 
and the ith lens to the jth lens respectively. The final Jaco- 
bian can be written in the form 



A : 



-fp21 



- 4>22 



(8) 



from which the components of the overall two-dimensional 
shear, 7, for weak lensing are 



71 = r 
and 

72 = y>2i 



' 1p2 



(9) 



(10) 



(In a weak shear field smoothed by the variable particle soft- 
ening and where the gravitational potential and its deriva- 
tives are well-behaved continuous functions, ip2i — ipi2-) The 
two-dimensional shear is 



7 = 71 + «72 



(11) 



(i = v 7 — 1), and the orientation of the major axis of the 
resulting elliptical image is 



-tan 



The effective convergence is 

« = ~C0u + ^22), 
and the final magnification is 
fj. = (det^) -1 = - 1 



' 7 



(12) 



(13) 



(14) 



Since the shear causes the axes to be stretched by fac- 
tors of 



a = (1 — k — I 7 |) _1 (major axis) 
and 

b — (1 — k+ I 7 |) _1 (minor axis), 



(15) 



(16) 



the imposed ellipticity can easily be calculated from the ele- 
ments of the Jacobian matrix. Observationally, the ellipticity 
of an image is frequently defined in terms of the tensor of 
second brightness moments, 



Qij 



J dflg f [J(fl)](fli -fli)(flj -Bj 

J de qi 



(17) 



where 1(9) is the surface brightness of the galaxy image at 
angular position 9, is the angular position of the centre 
of light, and qj[I(9)] is a weighting function in terms of the 
surface brightness. (See Blandford et al., 1991, for example.) 
Then one definition is the complex ellipticity 

Q11 - Q22 + 2iQia 



+ ' 



+ 2(QnQ 2 



(18) 



For elliptical isophotes, this definition is equivalent to 



1 - r 2i<t, 
1 + r 



(19) 



where r = b/a. 

The "reduced shear," g, evaluated for position 9 in the 
image, is defined by 



9(0) 



7(g) 
1 - k{0) 



(20) 



so that the transformation between the source and image 
ellipticities may be given by 



1 - ,9*e 



(21) 



for I 9 l< 1 (where the asterisk refers to the complex conju- 
gate). 

Then, in the case of weak lensing, for which k and 
I 7 |<C 1, I g |<C 1, so that, for low intrinsic source elliptici- 
ties, e ~ e (s ) + g. 

However, the intrinsic ellipticities of observed galaxies 
are generally unknown, so that the determination of the 
shear signal from individual images is impossible. For this 
reason, it is necessary to consider ensembles of galaxy images 
together, and to assume that the galaxies of each ensemble 
have random intrinsic ellipticities, so that the ensemble has 
zero net ellipticity. Whilst this is strictly not true, in high- 
redshift surveys in which the galaxies within each narrow 
cone may be widely separated, it serves as a good working 
approximation. A number of studies have been made into 
intrinsic correlations of galaxy shapes; Heavens, Refregier & 
Heymans (2000), for example, have shown that the intrin- 
sic correlation function for elliptical galaxies at z = 1 in an 
LCDM cosmology is only of order 10 -4 on angular scales of 
0.1' to 10', and is approximately an order of magnitude lower 
than the correlations expected from weak lensing. However, 
the intrinsic correlations are expected to exceed those arising 
from weak shear in shallow surveys, as found, for example 
by Brown et al. (2000), for sources with a median redshift 
of only 0.1. 

If then the ensemble of sources has zero net ellipticity, 

(e (s) ) = and 



(e) 



(22) 



where the u; are weight factors. Consequently then, in the 
case of weak lensing only, 



7 - 9 - (e) > 



(23) 



and the variances in both the shear and the reduced shear 
for a given angular scale are expected to be similar. 

The importance of the convergence for understanding 
the evolution of structure lies in its close association with 
the density contrast, <5(x), at position x, which is defined by 



*(*) = 



P(x) 



(24) 



By extension of the above equations for ipij (equation g) and 
k (equation |l3|), twice the value of the effective convergence 
in the direction 9 for a source at distance x a is 



2k(9,x s ) = 



D(x 3 )D(x s 



D(x s 



^1(V 2 



X7l g )<j>(9,x 3 )dx 3 
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D(x 3 )D(x s 



D(x s ) 



in which 

V 2 <t>(8,x 3 ) 



3Hi 5(0, x 3 ) 



a(x 3 ) 



'-^■V 2 4>{e,x 3 )dx 3 , (25) 



(26) 



W 2 [kD(x)<ff\ = 2Ji[kD(x)&\/[kD(x)'B], 



(33) 



where Ho is the Hubble parameter. Consequently, the ef- 
fective convergence represents a projection of the density 
contrast, and is proportional to the density parameter, fl m : 



2k(0,x s 



3Hi 



n 



D(x 3 )D(x s - x 3 ) 5(0, x 3 ] 



D(x s 



a(x 3 ) 



dx 3 .(27) 



It is also important to note how the two-point statistical 
properties of the shear and convergence are related. From 
the definitions of the individual shear tensors, U % (equa- 
tion 0) , the components 71 and 72 of the shear (equations ^ 
and and the effective convergence, k (equation the 
following expressions apply in Fourier space: 



7i (1) 
and 

72(1) 



1 2 



(28) 



(29) 



where li and I2 are the components of the wavevector 1, so 
that 



7i 2 (l) + 72 2 (l) = ft 2 (l). 



(30) 



Then it is clear that the power spectra for the shear, P-y(l), 
and the convergence, P K (l), are the same in the case of weak 
lensing. 

To obtain values for the shear (or convergence) vari- 
ances analytically, the convergence power spectrum is in- 
tegrated over all wavenumbers, using a filter function ap- 
propriate for the required angular scale, 9. Since the con- 
vergence is obtained from a projection of the density con- 
trast (equation ^7j) from the source redshift to the observer, 
the shear variance calculation requires a complete spatial 
and temporal description of the matter power spectrum, 
P$(k,x). This is a function of the real-space wavenumber, 
k (— l/(D(x)6), where D(x) is the angular diameter dis- 
tance for a radial distance x from the observer. Kaiser (1998) 
has determined general expressions for the angular power 
spectra of weak lensing distortions for different cosmolog- 
ical models, and has estimated the growth of this power 
with source redshift. As expected, the redshift dependence 
is much stronger in low density cosmologies and especially 
so in cosmologies dominated by a cosmological constant. 

Jain & Seljak (1997) give an expression equivalent to the 
following for the shear variances derived analytically from 
the matter power spectrum. 

/-CX3 

<7 2 )W = - x 367r 2 ft^ / fcdfc 
^ Jo 

x ( a- 2 (x)P s (k,x)g 2 (x)Wl[kD(x)-d]Ax, (31) 
Jo 

where x s is the radial distance to the source, 
D(x)D(x s - x) 



and 



D(x s ) 



(32) 



where Ji is the first Bessel function of the first order. The 
scale 1? is the angular radius of a circular window, so that the 
formula has to be transformed to express the shear variance 
on scales, 8, represented by square pixels. 

A complete description of the power spectrum on all 
scales (including the linear, quasi-linear and non-linear 
regimes) is necessary for the analytical approach, together 
with a detailed description of its evolution. In particular, 
non-linear effects on scales of order V may increase the 
amplitude of the convergence power spectrum by an or- 
der of magnitude. In addition, density fluctuations on scales 
smaller than about 10' contribute most strongly to the weak 
lensing signal, precisely where the non-linear evolution of 
the power spectrum is most in evidence. In the non-linear 
regime, the fitting formula? of Peacock & Dodds (1996), 
which extends the earlier work of Hamilton et al., (1991) 
for the evolution of the matter correlation function, may be 
used to map the non-linear wavenumbers onto equivalent lin- 
ear wavenumbers, and thus to evaluate the shear variance 
values. These fitting formulae include the stable-clustering 
hypothesis which assumes an invariant mean particle separa- 
tion on sufficiently small scales. Analytical programs based 
on this prescription are accurate to ~ 15%, depending on 
the cosmological model. 

For (circular) angular scales of # = 2' and 15', Jain 
& Seljak (1997) summarised their findings in approximate 
power-law expressions, which are equivalent to 

<7 2 )[tf = 2'(15')] oc 1 9- - 84 z s 1 ' 52 a 8 2 - 58(2 ' 00) ^ 20(1 - 36) (34) 

for LCDM cosmologies. The different indices for ag and fi m 
on the different scales enable the degeneracy between these 
parameters to be lifted when measurements are made in 
both regimes. The power of z s quoted is an intermediate 
value for the two angular scales. 

On the scales 1' < ■& < 30', the angular scale depen- 
dency was similar, and the powers of z s were close; the in- 
dex of z s decreased from 1.54 at 1? = 2' to 1.48 at 1? = 15'. 
As a result, Jain & Seljak (1997) found that, for LCDM 
cosmologies, 



<7 2 }(0) oc e-°- SA zY° 2 (l'<6»<30') 



(35) 



approximately. (1? and 9 are interchangeable in this context.) 

The numerically determined shear variance results 
quoted in this paper, and derived from the real-space values 
of the shear computed numerically, will be compared with 
this approximate expression. In addition, as a check of the 
validity of the numerical results, analytical values for the 
shear variance have also been computed directly for the pre- 
cise cosmology used in the numerical simulations. To achieve 
this the analytical program described in Van Waerbeke et 
al. (2001a) has been used, which is a quite general program 
for determining the weak lensing statistics in different cos- 
mologies. For the non-linear evolution of the power spectrum 
(and the determination of variances), the fitting formula? of 
Peacock & Dodds (1996) are used. For higher-order statis- 
tics, such as the weak lensing skewness (not computed in 
this work), the code also computes the evolution of the bis- 
pectrum at all scales, based on fitting formula? derived by 
Scoccimarro & Couchman (2001) in numerical simulations. 
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3 NUMERICAL PROCEDURE 

Couchman et al. (1999) describe in detail the algorithm for 
the computation of the elements of the matrix of second 
derivatives of the gravitational potential in cosmological TV- 
body simulations. It computes all of the six independent 
component values of the three-dimensional shear at each of 
the selected evaluation positions. The rms errors in the com- 
puted shear component values are typically ~ 0.3%. 

The algorithm uses a variable particle softening which 
distributes the mass of each particle throughout a radius 
which depends on its specific environment. The actual values 
of the softening parameters are precisely as described in Bar- 
ber et al. (2000), with the minimum value set to 0.1/i _1 Mpc 
throughout, where h is the Hubble parameter expressed in 
units of 100km s _1 Mpc _1 . 

In the computation of the shear, the code uses the pe- 
culiar gravitational potential, cj>, through the subtraction of 
a term depending upon the mean density. This ensures that 
only light ray deflections arising from departures from ho- 
mogeneity apply, and is equivalent to requiring that the net 
total mass in the system be set to zero. The algorithm auto- 
matically includes the contributions of the periodic images of 
the fundamental volume in computing the three-dimensional 
shear at any location, thereby essentially creating a realisa- 
tion extending to infinity. 

Since the algorithm works within three-dimensional 
simulation volumes, rather than on planar projections of the 
particle distributions, angular diameter distances to every 
evaluation position can be applied. In this work it has been 
assumed that the angular diameter distance varies linearly 
through the depth of each simulation volume. 

The code has been applied to the cosmological TV-body 
simulations of the Hydra Consortium]!] produced using the 
'Hydra' TV-body hydrodynamics code (Couchman, Thomas 
& Pearce, 1995). Simulations of the LCDM Dark Matter 
only cosmology were used with Q m = 0.3, Ao = 0.7, power 
spectrum shape parameter V = 0.25 and normalisation, as, 
on scales of 8/i _1 Mpc of 1.22. The number of particles, each 
of mass 1.29 x lO 11 /)." 1 solar masses, was 86 3 and the simu- 
lation boxes had comoving side dimensions of 100/i _1 Mpc. 
The simulation output times were chosen so that consecutive 
simulation boxes could be abutted. A total of 48 boxes to a 
redshift of 3.57 in the LCDM cosmology were used. To avoid 
obvious structure correlations between adjacent boxes, each 
was arbitrarily translated, rotated (by multiples of 90°) and 
reflected about each coordinate axis, and in addition, each 
complete run was performed 20 times, so that averages of 
the final statistics were determined to represent the required 
results. 

To follow the behaviour of light rays from distant 
sources through the simulation boxes, and obtain distribu- 
tions of the properties at z = 0, a set of light paths was 
constructed emanating from the centre of the front face of 
the 2 = box and ending in a regular square array of loca- 
tions at the plane of the chosen source redshift. 

A total of 14 source redshifts were selected to give good 
statistical coverage of the redshifts of interest. These were 
redshifts of z s = 0.10, 0.21, 0.29, 0.41, 0.49, 0.58, 0.72, 0.82, 
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0.88, 0.99, 1.53, 1.97, 3.07 and 3.57. They corresponded to 
the simulation box redshifts and were chosen to be close to 
redshifts of 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.5, 
2.0, 3.0 and 3.5. In this paper the source redshifts will be 
referred to loosely as the latter approximate values, although 
in the determination of redshift dependences, etc., the actual 
redshift values were used. 

Using a total of 317 x 317 lines of sight, the angular size 
of the minimum particle softening is comparable to or less 
than the angular separation of the adjacent lines of sight 
(0'.49) for all redshifts greater than 0.14. Consequently, to 
account for the larger angular size of the minimum softening 
at low redshifts, and also for the allowed range of particle 
softening scales above the minimum value, a resolution limit 
of 2' has been adopted for the data analyses. In addition, the 
angular size of the gravitational force softening used in the 
generation of the simulation boxes is below the line-of-sight 
separation for all redshifts greater than 0.05. For a source at 
redshift 1 in the LCDM cosmology, maximal gravitational 
lensing occurs for a lens at redshift 0.36, and the angular sep- 
aration of adjacent lines of sight is approximately the same 
value as the angular size of the minimum particle softening 
at that redshift (0'.47). 

The total field-of-view for the set of lines of sight was 
2.6° x 2.6°, and this completely fills the near face of the 
simulation box at redshift 1.0. 

To establish the locations for the evaluation of the shear 
on each of the lines of sight, first a regularly spaced (coarse) 
set of 50 locations was laid down on each line in each simu- 
lation box. Then additional locations were computed at po- 
sitions where the gravitational potential was changing most 
rapidly, so that the potential field could be well-sampled. 
To establish these locations, the particles were assigned to 
volumes determined by a 10 x 10 grid within each box. From 
each of the coarse line of sight locations, the separations to 
the particles within the local grid volume and the nearest 
neighbouring grid volumes was determined. If a separation 
was less than the line-of-sight separation, a new evaluation 
location was established on the line-of-sight, with coordi- 
nates corresponding to the particle's position in the xz di- 
rection (radially from the observer). All the evaluation lo- 
cations (coarse and new) along each line in each box were 
then sorted, labelled and counted so that the programme to 
integrate the values along the lines of sight would operate 
in the correct order for the correct number of locations. 

Following the shear computations at all of the locations 
on all the lines of sight in all the simulation volumes, the 
second derivatives of the two-dimensional effective lensing 
potentials were obtained from the three-dimensional shear 
values by integration, in accordance with equation 0. The 
integration was made in step-sizes determined by the sepa- 
ration of adjacent evaluation locations on each of the lines 
of sight, and so was different for every pair of points. The 
integration algorithm was set to run from each of the cho- 
sen source redshift planes, along each of the lines of sight to 
z — 0, and values for the elements of the shear matrix, and 
thus the Jacobian matrix, at the observer were obtained for 
each of the lines of sight. From these data, all the required 
weak lensing statistics were obtained for each line of sight 
and for each source redshift. 

The full procedure, from the computation of the three- 
dimensional shear values at all the evaluation locations to 
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finally obtaining the Jacobian matrices at z = 0, involved 
precisely the same approximations as described fully in Bar- 
ber et al. (2000). In addition, it should be noted that, be- 
cause the lines of sight project radially from the observer at 
z = 0, some lines of sight pass outside the confines of the 
simulation volume beyond a redshift of 1. In these cases, the 
periodicity of the boxes has been used to reposition the lines 
of sight at equivalent locations within the volumes. The area 
of the far face of the most distant simulation box (at red- 
shift 3.6) is 1.3° x 1.3°. Consequently, the procedure is not 
expected to introduce significant effects on the computed 
variances on the scales of interest here, which are up to a 
maximum of 32'. 



4 RESULTS 

The individual real-space values for the shear, 7, and the 
reduced shear, g, were computed from the final Jacobian 
matrices obtained in each run for each line-of-sight and for 
each source redshift. These data sets were then separately 
convolved with a top-hat smoothing function of the required 
scale-size, and the statistical variance values obtained on 
those specified scales. The scale-sizes chosen for the top-hat 
smoothing were 2' (consistent with the resolution limit of the 
numerical procedure), 4', 8', 16' and 32'. The computed 
values from each of the N runs were then averaged, and the 
errors on the means of Xa/yN determined. 

Figure |l| shows the shear variances computed in this 
way (without the error bars for clarity) for all the source 
redshifts, with the exception of z s = 0.1, which is too low to 
be seen clearly. The figure clearly emphasises the redshift de- 
pendence of the results and suggests that a good knowledge 
of the redshift distribution of the sources observed in surveys 
is essential to interpreting the shear signal correctly. The rel- 
ative closeness of adjacent curves separated by intervals of 
only 0.1 in redshift also suggests that the shear signal from 
sources with a small spread in redshift may be adequately 
described by the shear signal expected from sources at their 
median redshift. 

However, it is not clear that the shear resulting from 
sources with a significant redshift distribution will be repre- 
sentative of the shear from sources at their median redshift. 

A number of recent measurements of cosmic shear have 
been made, as mentioned in the Introduction. Following 
these measurements, Kaiser et al. (2000) and Bartelmann 
& Schneider (2001) have plotted them on a single diagram 
and compared the results with the predicted values for the 
shear variances as suggested by Jain & Seljak (1997), whose 
work is described above. As can be seen from the diagrams 
in these references, the cosmic shear signal resulting from all 
the observed measurements (from different telescopes, filters 
and cameras, and different fields of view and data analysis 
techniques), appears to lie very close to the formerly pre- 
dicted values for sources at a redshift of 1. 

The numerical values computed here for source redshifts 
of 0.4, 0.6, 0.8, 1.0 and 1.5 are now plotted in figure ^| 
together with the analytical values determined using Van 
Waerbeke et al.s (2001a) code, described above; also shown 
are Jain & Seljak's (1997) predicted values (transformed to 
square pixel areas) for source redshifts of 1, using their ap- 
proximate expression for general LCDM cosmologies. The 
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Figure 1. (-y 2 (6)) for source redshifts 0.2 (lowest curve), 0.3, 
0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.5, 2.0, 3.0 and 3.5 (uppermost 
curve) . 

advantage of Van Waerbeke et al's (2001a) prescription is 
that the precise cosmological parameters and source red- 
shifts used in the numerical work have been used. 

The numerically computed values for z 3 — 1 are remark- 
ably close to the values predicted analytically by both Jain 
& Seljak (1997) and Van Waerbeke et al.s (2001a) prescrip- 
tion. In particular, the analytical results computed here give 
strong support to the numerical results, and the agreement is 
particularly good at low redshift and at intermediate angu- 
lar scales. The largest discrepancies occur only for sources 
at high redshift (beyond about 1.5) and at angular scales 
comparable with and below the resolution limit of the nu- 
merical data, and where breakdown of both the numerical 
and analytical procedures may also be expected. The values 
from Jain & Seljak's (1997) approximate general expression 
differ somewhat from both the numerical and analytical re- 
sults reported here for source redshifts other than 1, so that 
it appears qualitatively that the results here would show 
a different redshift dependence than the former analytical 
predictions. 

The redshift dependence of the shear variance on the 
chosen angular scales is now plotted in figure ^. 

Since the redshift relationship for the shear variance 
and its scale dependence are of fundamental importance to 
observational measurements of the cosmic shear signal, the 
numerical results from figures [j], ^ and |^ have been fitted 
mathematically. It is assumed that the shear variance can 
be expressed in the simple form 

tf){e,z a )=a{6)z b 3 w , (36) 

and indeed, for the redshift range z s < 1.6, this form de- 
scribes the shear variances well. 

Over the range of angular scales 2' to 32', the coefficient 
a(8) can be expressed as 

a(6) = (1.05 ± 0.05) x io-^" 1 ' 12 * 03 , (37) 

for 8 in arcminutes. At the low end of the range of scales, the 
index b(6) falls very slightly as 6 increases, and then remains 
almost constant to beyond 32'. Hence, for the entire range 
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Figure 2. The numerical shear variances and analytical values 
computed here, based on an identical cosmology and the same 
source redshifts. Source redshifts of 0.4 (lowest pair of curves), 
0.6, 0.8, 1.0 and 1.5 (uppermost pair of curves) are plotted, re- 
spectively. For comparison, the analytical predictions of Jain & 
Seljak (1997) for sources at a redshift of 1 are also plotted, based 
on their approximate general expression for LCDM cosmologies. 



Figure 4. (7 2 (0)} vs. z s on logarithmic axes, to show the gradu- 
ally declining slopes at high redshift for the measured shear vari- 
ances. The curves are for angular scales of 2' (uppermost curve), 
4', 8', 16' and 32' (lowest curve). 
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Figure 3. (7 2 (#, z a )) as a function of source redshift for angular 
scales of 2', 4', 8', 16' and 32'. 



of scales 2' to 32', b(8) can be described well by the constant 
value 

b{9) = 2.07 ±0.04. (38) 

The relationships amongst {■y 2 }, and z s may be seen 
more clearly in figure ^ for the numerical data, which has 
logarithmic axes. At low redshifts the slopes of the curves 
have a very small 6 dependence, and, for redshifts greater 
than about 1, the curves become clearly less steep, indicating 
that the redshift dependence falls well below 2 at high source 
redshifts. 

In Section 2, it was shown that, in the case of weak 
lensing, both the shear, 7, and the reduced shear, g, will 
be approximately equal, so that the variances in both these 



.2 o 





1 1 1 1 1 1 
shear 




- \ 






4. \ 


reduced shear 




- \ \ 












\ x 






\\ 






\\ 






\\ 












V" 














z =2.0 





























5 10 20 

Angular scale (arcminutes) 



Figure 5. The variances in the reduced shear and the shear for 
comparison, for source redshifts of 0.5, 1.0 and 2.0. 



quantities will also be similar. Indeed, this approximation is 
used observationally to estimate the values of the shear from 
the observed ellipticities. However, the equality holds only 
in the weak lensing limit. The numerical method used here 
for computing the weak lensing statistics enables a direct 
comparison to be made between the variances in 7 and g. 

Figure H compares {g 2 }(0) with (7 2 }(0) for redshifts of 
0.5, 1 and 2. For low source redshifts, the curves are almost 
identical (except for the smallest of angular scales). This is 
as expected, because this is the regime of the weak lensing 
limit. However, departures between the two quantities be- 
come increasingly obvious as the source redshift is increased. 

The redshift dependence for {g 2 } on the different angu- 
lar scales is plotted in figure ^[ 

The functional form of (g 2 )(8,z a ) can be written in a 
similar way to that for (■y 2 }(8, z s ). By analogy with equa- 
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Figure 6. (g 2 (0,z s )) as a function of redshift for the angular 
scales of 2', 4', 8', 16' and 32'. 

tion the variance in the reduced shear can be written 

as 

(g 2 ){9,z a )=c{9)z d s W (39) 

in the range z s < 1.6, for 2' < 9 < 32'. 

The coefficient c(9) can be expressed as 

c{9) = (1.2 ± 0.1) x io- 3 ^- 1 - 21 * 03 . (40) 

As for b{9) in the expression for the shear variance, d{9) 
here falls very slightly as 9 increases at small angular scales, 
but then remains almost constant throughout the range of 
angular scales. Explicitly, for 2' < 9 < 32', d(9) can be 
expressed as the constant value 

d(9) = 2.01 ± 0.04. (41) 

The variance in the reduced shear for the entire redshift 
range is plotted with logarithmic axes in figure |^ to indicate 
both the 9 and z s dependencies of the results. 

As found for the shear variance, figure 0for the reduced 
shear, shows how the slopes of the curves at the low redshift 
end have a slight 9 dependence. However, unlike the shear, 
at high redshifts the slopes become more noisy and gener- 
ally steeper, even though the slopes do decline very gently 
throughout the redshift range up to about 2. This is most 
likely a result of the sensitivity to the convergence, apparent 
from equation which becomes important either on small 
angular scales or at high redshifts, where the lensing effects 
are strongest. 



5 DISCUSSION AND CONCLUSIONS 

The main results established in this paper are the depen- 
dences of the shear variance and the variance in the reduced 
shear on source redshift and angular scale. The numerical 
data obtained has been validated by analytical computa- 
tions using Van Waerbeke et al.'s (2001a) program applied 
to an identical cosmology and the same source redshifts. 

For the shear variance, for source redshifts up to 1.6, 
one can write, without the error bars for clarity, 



Figure 7. log({g 2 (0))) vs. log(,z a ), to show the gradually declin- 
ing slopes throughout most of the redshift range for the reduced 
shear variances on angular scales of 2' (uppermost curve) , 4' , 8' , 
16' and 32' (lowest curve). 

{-y 2 )(9,z a ) = 1.05 x lO^lT 1 - 12 ^ 07 (42) 

for 2' < 9 < 32'. I.e., the redshift dependence is close to 
z 2 , and stronger than the previous analytical predictions 
of Jain & Seljak (1997) who found a dependence close to 
zl' 52 . Whilst the above equation fits the numerical data well 
within the range specified, the redshift dependence clearly 
begins to decline for sources beyond a redshift of 1, so that 
the index of z s falls substantially below 2 at high redshift. 

For the variance in the reduced shear, again it has been 
found that the redshift dependence is close to z 2 , and for 
z a < 1.6 it can be described by 

(g 2 )(9,z s ) = 1.2 x Mr 8 *" 1 - 21 **- 01 (43) 

for 2' < 9 < 32', again without the error bars for clarity. At 
higher redshifts, although there is an underlying tendency 
for the curves to become less steep, the reduced shear is 
noisy and sensitive to values of the convergence, so that the 
curves steepen again at the high redshift end. 

It was mentioned in the Results section that Kaiser et 
al. (2000) and Bartelmann & Schneider (2001) have plotted 
recent observational determinations of the shear variances 
on a single diagram and compared the results with the pre- 
dicted analytical values of Jain & Seljak (1997). These plots 
show the observational results to lie close to the predicted 
curve for sources at a redshift of 1. It should be noted at this 
point that the values determined numerically here, and also 
analytically using Van Waerbeke et al.'s (2001a) program, 
concur well with the analytical values plotted for source red- 
shifts of 1. However, the observational results plotted have 
galaxies with distributions in redshift; for example, Bacon et 
al. (2000) quote a median redshift of z s = 0.8 ± 0.2 for their 
sample, Van Waerbeke et al. (2000a) refer to a peak redshift 
of 0.9 for their galaxies, Kaiser et al. (2000) has an "effec- 
tive" redshift of 1.0, and the survey of Maoli et al. (2001) 
has a broad redshift distribution which peaks at a redshift 
of 0.8. 

Whilst Jain et al. (2000) claim that a distribution of 
sources with a mean redshift of 1 gives rise to an amplitude 
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for the shear variance only 10% different from the case where 
all the galaxies are assumed to lie at z s = 1, the contention 
here is that it is necessary to have a clear understanding of 
how the shear signal relates to the redshift of the sources, 
and how a distribution of source redshifts may further influ- 
ence the conclusions. As an example, there are differences of 
~ 10% - 20% between the shear variance values for sources 
separated by a redshift interval of only 0.1 on angular scales 
of 2'. Attempts to constrain the cosmology precisely will 
founder if the redshift dependence reported here is not taken 
into account. According to equation ^ an uncertainty of 
about 0.025 in the density parameter also results if there is 
a 10% uncertainty in the shear variance. In addition, there 
are clear differences between the variances in the shear and 
the more directly measured reduced shear, particularly at 
high redshift and on small angular scales where the weak- 
lensing regime may not be applicable. These factors point 
clearly to the need for studies into the strength of the cosmic 
shear signal for sources distributed in redshift, and will be 
investigated in a future paper. 

The faster redshift dependence reported here (~ z 2 ' 07 
compared with ~ zl' 52 previously predicted) is not easily ex- 
plained, as there are assumptions in both approaches. In the 
numerical method one could point to the discontinuities in 
structure as one passes from box to box, although the effect 
of this is reduced statistically by employing a large number 
of runs (20 in this work). The particle softening and resolu- 
tion limitations might also be a factor in the simulations. In 
the analytical approach, the main difficulty arises in trying 
to describe accurately the evolution of the power spectrum 
in the non-linear regime. Here mapping techniques using fit- 
ting formulae, accurate to ~ 15%, are used to relate the 
linear and non-linear spectra, and assumptions such as the 
stable-clustering hypothesis are used. Also, Jain & Seljak's 
(1997) published predictions are the result of summarising 
data from different cosmological models, and then fitting the 
data to approximate formulas to express the power-law de- 
pendencies of the statistics. Improvements have been made 
more recently however in analytical procedures, for example, 
Seljak (2000) and Peacock & Smith (2000) have developed 
a model for the non-linear evolution of the power spectrum 
based on the random distribution of dark matter haloes, 
modulated by the large-scale matter distribution. 

The main advantage of the numerical method is that 
the shear values and variances are determined directly from 
the real-space data, rather than from integration of a power 
spectrum which is less well understood on the scales of inter- 
est. It is also expected that numerical results may be used in 
the future to test the validity of developments in analytical 
procedures, which will be increasingly needed in view of the 
expected accuracies in the next generation of cosmic shear 
measurements. For the interpretation of these observations, 
fast analytical procedures will be more practical, since it 
will be unrealistic to generate large numbers of numerical 
simulations to cover a realistic parameter space. 

A key influence on the redshift dependence of the shear 
variance lies in the functional form of Q (x), the angular di- 
ameter distance function. Q(x) occurs in equation |3l| where 
it multiplies the matter power spectrum to produce the effec- 
tive convergence power spectrum. In the numerical method 
here it also plays a crucial role in establishing the elements of 
the final Jacobian matrices, from which the convergence and 



shear values are obtained. In the numerical work here, the 
"filled beam" values for Q(x) have been used as it is assumed 
that the light is always passing through a region of smoothed 
density. (See Barber et al., 2000, for a full discussion of an- 
gular diameter distances in inhomogeneous universes.) The 
form of G(x) at high redshifts may also help to explain the 
reducing index of z 3 , because the redshift value at the peak 
of the Q(x) curve rises less quickly as the source redshift is 
increased. 

In this paper, no attempt has been made to simulate 
the signal which might arise from typical noisy data. The 
intrinsic ellipticities of the background galaxies is a source 
of random noise, and in addition, measurement errors of very 
faint galaxies introduce errors in the ellipticity estimates. 

The intrinsic ellipticity correlations of galaxies in close 
proximity to each other produce a real signal which has to 
be disentangled from the cosmic shear correlations from the 
large-scale structure along the line-of-sight. Heavens et al. 
(2000) have attempted to determine the intrinsic correla- 
tions by modelling the shapes and angular momentum vec- 
tors of galaxies in a way which reflects the shape and an- 
gular momenta of their dark matter halos in iV-body sim- 
ulations. They find that the intrinsic ellipticity correlation 
function for elliptical galaxies is an order of magnitude be- 
low the expected weak lensing signal on scales of 0.1' to 10' 
for sources at z s — 1 in a LCDM cosmology. Brown et al. 
(2000) have analysed real data from the SuperCOSMOS Sky 
Survey for galaxies with a low median redshift of z s = 0.1. 
They find that the ellipticity variance is approximately two 
orders of magnitude higher than the expected weak lensing 
signal throughout the range 1' to 100'. As the source red- 
shifts are low, they claim that their result represents the real 
intrinsic correlations since the weak lensing signal would be 
expected to be small. 

Taking the above two reports into account, there will be 
some intermediate redshift at which the weak lensing signal 
overtakes the intrinsic one. Consequently, whilst the result 
reported here presents the pure redshift dependency of the 
weak lensing signal, it is clearly advantageous to study good 
data at high redshift in order to measure the uncontami- 
nated cosmic signal. 

This represents a further argument in favour of hav- 
ing good redshift information in a galaxy survey designed 
to measure cosmic shear. It is anticipated that forthcom- 
ing deep weak lensing surveys, such as that proposed for 
the VISTA telescope, will also provide detailed photomet- 
ric redshift information to enable the surveyed galaxies to 
be binned in redshift intervals. In this way it is anticipated 
that the real cosmic shear signal may be interpreted more 
correctly. 
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